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Abstract 

In this paper we present a Mathematica notebook for computing 
nonsymmetric and interpolation Macdonald polynomials. We present 
the new recursive generation algorithm employed within the notebook 
and the theory required for its development. We detail the contents 
of the notebook and conclude with a couple of applications of the 
notebook. 

1 Introduction 

The nonsymmetric Macdonald polynomial E v (z; q,t) and its generalisation - 
the interpolation Macdonald polynomial E*(z; q, t) - have found applications 
in mathematical physics, combinatorics and representation theory [12J, [TIJ El 
E3]. 

A feature of the aforementioned polynomials is that they allow for explicit 
computation. Having this knowledge provides an opportunity to experimen- 
tally seek new properties and to check analytic work. Within this paper a 
computer software program that was developed for these purposes is detailed. 

The nonsymmetric and interpolation Macdonald polynomials are gener- 
ated via recursive generation formulas. This paper begins with the theory 
required to introduce these formulas (Section [2]) and then later shows how 
these formulas can be used to construct a recursive generation algorithm for 
the generation of the nonsymmetric and interpolation polynomials (Proposi- 
tion |2]). Section E3] details the functions within the Mathematica notebook 
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containing this algorithm and the paper concludes with two possible appli- 
cations of the notebook. 

2 Foundation Theory 

2.1 Compositions 

The polynomials discussed in this paper are labelled by compositions. We 
define a composition to be an n-tuple 77 := (771, . . . ,i] n ) of non-negative in- 
tegers. Each rji is called a component, the length of a composition is the 
number of components it contains and the sum of the components is called 
the modulus and denoted by \rj\. Each composition 77 corresponds to a unique 
partition 77+ obtained by rearranging the components of 77 so that they are 
nondecreasing. 

Two important operators that act on compositions are the switching and 
raising operators. We have the switching operator Sj acting on compositions 
according to 

Sifa, ...,1Ji, Tfi+i, ...,7ln):= (t]i, . . . , 1Ji+i, T)i, . . . , Ifa), % = 1, . . . , n - 1 

and the raising operator $ which has the action 

$77 := (T) 2 ,...,r} n ,Tii + 1). (1) 

It can be seen that every composition r) can be recursively generated from 
the all zero composition (0, . . . , 0). In Section 1X21 we construct an algorithm 
to recursively generate any composition from 77 from (0, . . . , 0) using the 
least number of operators. It is this algorithm that allows us to obtain our 
main result, the recursive generation algorithm for the nonsymmetric and 
interpolation polynomials (Proposition [2]). 

We now proceed to the required polynomial theory. 

2.2 Polynomial theory 

The nonsymmetric Macdonald polynomial [7J and interpolation Macdonald 
polynomials [H] are most commonly defined by their eigenfunction and van- 
ishing properties, respectively. Here we take an alternative approach and 
define them via their respective recursive generation operators. 
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In the previous section we noted that every composition r] can be recur- 
sively generated from (0, . . . , 0) using the operators and In a similar 
way the nonsymmetric Macdonald polynomial E v (z;q,t) := E v (z) and the 
interpolation Macdonald polynomial E*(z;q,t) := E*(z) can be recursively 
generated from £ , ( 0i ... i o)(-2; Q, t) = 1 and Ef -,(z;q,t) = 1, respectively, us- 
ing two elementary operators. We begin by introducing the switching-type 
operators for the two polynomials; namely the Demazure-Lustig operator 7] 
for the nonsymmetric polynomials and the Hecke operator Hi for the inter- 
polation polynomials. These operators relate E v (z) and E SiV (z), and E*(z) 
and E*. n (z), respectively. 

Each switching-type operator is a realisation of the type- A Hecke algebra, 
an associative unitaal algebra over Q(t) generated by elements hi, ... , 
and subject to the relations 

hihi + \hi = hi + ihihi + i, for 1 < i < n — 2 

hihj = hjhi, \i- j\> 1 (2) 
(hi + l)(hi-t) = 0. 

The switching type operators are defined by 

Ti ;=t+ tZi ~ Zi+1 { Si -l) and ff,:=t+ ^~ tem (s,-l), (3) 

where here, the switching operator Si is defined to act on functions according 
to 

\Sif)\Zlj • • • ) %ii Zi+li • • • i %n) ■ /(^lj • • • i ^i+li %ii • • • i ^n) ■ 

These operators act on their respective polynomials according to [T8] 
TiE v {z) = { tE v (z), rji = r] i+1 

Z>r)\Z) H Q r. ^Siri{Z), Tji > Tli + i, 



and (SHE] 



(4) 



t-i 



1-^(9.*) 



(5) 



i_ { -i (9 ,i)^l z ) + tosf ^vA^ Vi>Vi+i- 
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In flU and © S itV (q,t) := r)Jr) iJrl , with 
where 

W : = #{j < ^ > + #{j > > Vi}- 

The raising type-operators $ g and $* of the nonsymmetric and interpolation 
polynomials, respectively, are given by pQ 

$ g := z n T-\ . . . Tf 1 = t^ n T n ^ . . . T lZl Tr_\ . . . Tf 1 (6) 

and 

$* := (z n - t~ n+1 )A. (7) 

In ([6]) the operator is related to T, by the quadratic relation in fl2]) and 
given explicitly by 

rr 1 ■= r l - 1 + r 1 ^, 

and in ^ A/(zi, . . . , z n ) := f{z n /q, z u . . . , 

The raising type operators act on the Macdonald polynomials according 

to DP 

%E,(z)=t-#l i>1 ^} E ^ z) 

and pi] 

In Section 13.31 we provide details of the algorithm that implements these 
formulas to recursively generate any nonsymmetric or interpolation Macdon- 
ald polynomial. 

Below are some examples on nonsymmetric and interpolation polynomi- 
als, for further details on the polynomials we refer the reader to [I]. 



fT, - -y3 , t-l 2 1 , (g+lXt-l) 1 2 
^(0,3) \Z) - z 2 + jHZi z l z 2 H q H-\ Z l Z 2 

&(2,1)\Z) — Z X Z 2 + -^ zr Z l Z 2 

E(i,2)(z) = z\z\ 

E (1,0)( Z ) = z l + ^Zl z 2 

E*a,i)( z ) = z 2~\- 



- t Z 2 + W 

_ <?* 2 -l 

t{qt-X) 
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Before considering the recursive generation algorithm we discuss the lim- 
iting properties of the nonsymmetric and interpolation Macdonald polyno- 
mials. 

2.3 Limiting properties 

The following diagram indicates that the nonsymmetric Macdonald polyno- 
mials can be obtained from the interpolation Macdonald polynomials through 
a process of homogenisation, that is taking only the terms of top degree from 
the interpolation polynomial extracts the nonsymmetric polynomial. Fur- 
thermore, the diagram indicates that there are many additional families of 
polynomials that can be obtained from the nonsymmetric and interpolation 
Macdonald polynomials. It is this feature that allows the computational work 
discussed in this paper to immediately yield algorithms for the generation of 
these other families of polynomials. 

Through symmetrisation the nonsymmetric Macdonald polynomial E v (z; q, t) 
can be reduced to symmetric Macdonald polynomial P v +(z;q,t). Addition- 
ally, the symmetric, nonsymmetric and interpolation Macdonald polynomials 
reduce to the corresponding Jack polynomials, P v +(z; a), E v (z; a), E*(z; a), 
in the limit t = q 1 ^, q — > 1, which themselves reduce to the Schur and 
zonal polynomials, s v +(z), Z n +(z), by setting a to 1 and 2, respectively. 
Lastly, we have that the symmetric Macdonald polynomials reduce to the 
Hall-Littlewood polynomials P v +(z;t) in the limit q — 0. 
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S v + (z) Z v + (Z) 



The definitive reference for the symmetric Macdonald polynomials is [15J. 
This book also contains details on the Schur and Hall-Littlewood polynomi- 
als. An accessible reference for the Jack polynomials and zonal polynomials 
are [TU] and [TS], respectively. For details on how the symmetric Macdonald 
polynomials can be obtained from the nonsymmetric polynomials we refer 
the reader to [T7] . 



3 Computer Programming 
3.1 A brief history 

Computer software that generates the polynomial families listed above is 
quite limited. One of the few programs available is designed to compute 
symmetric Jack polynomials [9]. The software, written in Maple, uses recur- 
rence formulas to evaluate the coefficients K' (a) G Q(«) in the expansion 

P v +(z;a)= K \+ v +( a ) m \+( z )i 

\+<ri+ 

where the m\+(z) are symmetric monomial functions (see for example [15]). 
The study in [9] discusses the implementation of the software and the run 
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times of particular functions. They also provide details for the generation of 
the generalised classical Hermite, Laguerre and Jacobi polynomials. 

A further study in the generation of symmetric Jack polynomials is by 
Demmel and Koev |8] . They use the expansion formula for a Jack polynomial 
in one of its variables to obtain a more efficient evaluation for sets of Jack 
polynomials than was known previously. 

Although there are known methods for generating the nonsymmetric and 
interpolation Macdonald polynomials, for example the Rodrigues formulas 
[20] . software that generates nonsymmetric and interpolation polynomials 
appears to be nonexistent in the literature. It is our aim to initiate momen- 
tum in this area. 

We note that the methods presented here would be dramatically im- 
proved in efficiency if the algorithms introduced by Demmel and Koev were 
generalisable to nonsymmetric theory. Although progress has been made 
using computer-generated coefficients, confirming the ability to generalise 
the algorithms, the formulas cannot be run efficiently until the nonsymmet- 
ric equivalent of equation (6.2) of [8] - - a type of dual Pieri formula — is 
known. 

We now present the algorithms implemented in our Mathematica note- 
book. 

3.2 Recursively generating compositions 

As stated earlier every composition 77 can be recursively generated from 
the all zero composition (0, . . . , 0) using a (non-unique) sequence of raising 
and switching operators. For example, rj = (0, 2, 1) can be generated from 
(0,0,0) by SiS 2 SiS2$S2$s 2 $(0, 0, 0) or more efficiently by s 2 $Si$$(0, 0, 0). 
In polynomial classes for which switching and raising operators exist analo- 
gous methods can be used to generate polynomials. Consequently once we 
construct an algorithm that generates any composition recursively we auto- 
matically obtain an algorithm for the polynomials. 

Our aim is to construct an algorithm to recursively generate any compo- 
sition 77 from (0, . . . , 0) using the least number of operators. We first observe 
that the raising operator $ must be used \rj\ times. Since the raising operator 
acts on a composition by increasing the value of the component in the first 
position by one, appending it to the end of the composition and shifting each 
other component back one position, to minimise the number of operators we 
must always increase the value of the leftmost component requiring raising. 
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A systematic way of doing this is to apply the raising operator to build all 
components greater or equal to a specific size only using the switchings to 
move the leftmost component needing raising to the first position. We note 
that this method is quite similar to the Rodrigues formulas construction [20J. 
Using this method we naturally construct the composition (r] + ) R , where 

V R ■= (Vn,Vn-l, ■ ■ .,7fr). 

Due to the nature of $ there is no possible way to construct a composition 
containing each component of rj using fewer operators. 

To reorder (r] + ) R minimally we switch each component into its correct 
position beginning with either rji, 772, . . . or rj n , rj n -i, ■ ■ ■■ We choose to start 
with repositioning rj n . By always choosing the closest component of the 
unordered composition we ensure that the number of switches is minimal. 

Proposition 1 Define l ni := ^{rjj < i}, g Vj i := i^{f]j > 2} and 

^Wj -= i}r),i—li •■■)!; lri,i—l + 1, • • • , 2, . . . , Z^j—l + g^j 1, . . . , Qr],i)- 

Define 



where s^...^) := s it . . . Sj x . Define 



i-1 

p Vti := max j j < i : JJ s a , {ruk) (r] + )f = rji 

k=l 

and 

s (p v ,i,...,i-l)> Prj,i < 1 



Srr' 



1 ; Vr\,i i- 



The minimal length sequence of operators that transforms (0, . . . , 0) to tj is 
That is 

s <t'(77,2) • • • V(i),n)'"r ) ,max(»)) • • • r »?,l(0> • • • , 0) = 7?. (8) 
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Proof We prove (jSD in two steps. We first show by induction that 

r^axM-.-rvtO,...^) = (T)+) R . (9) 

By the definition of r Vt i it is clear that r v ^(0, . . . , 0) produces a composi- 
tion of the form (0, . . . , 0, 1, . . . , 1) where the number of l's is equal to the 
number of components of rj that are greater or equal to 1. Suppose before 
applying r Vt k+i we have generated the correct number of components with 
value 0, 1, . . . , k — 1, that is we have constructed a composition of the form 
((i] + )f , ,...,(ri + )^,k...,k) where the number of fc's is equal to the num- 
ber of components of rj that are greater or equal to k. Quite obviously 
r„, fc+1 ((r/+)f , . . . , ( V + )f, k...,k) = (fa+)f , . . . , (7j+)f , k...,k, k+1, . . . , fc+1) 
where the number of fc's equals the number of components of r] equal to k 
and the number of (k + l)'s equals the number of components greater than 
or equal to k + 1. By induction holds. The final task is to show that 

S*'( V ,2) ■ ■ ■ S a '( V ,n)(r] + ) R = V- (10) 

This result follows immediately from the definition of as quite clearly 
a'(rj,i) successively permutes each rji into the correct position. The fact that 
the total sequence of operators is of minimal length follows from the action 
of $ and the inability to generate a composition with components 771, . . . , r\ n 
more economically than what is specified by (Q, the permutation that places 
each component into its correct position can not be improved either. □ 

We note that further evidence showing the permutation in ffTU]) is minimal 
is the comparison of its length to the minimal permutation lu^uT+^h, where 
u v is the minimal length permutation such that o;~ 1 (r/) = r/ + . Due to the 
different structures of the permutations we use the computational evidence 
to support our claim. 

To provide additional clarity to the algorithm we show how (4, 1, 2, 1) is 
generated using the above operators. 

Example 1 We construct the composition (4, 1, 2, 1) recursively from (0, 0, 0, 0). 
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We first construct (1, 1, 2, 4) using the operators r Vti . 

r V A r v,3 r v,2r v ,i(0, 0, 0, 0) = r V)4 r V}3 r V} 2^(0, 0, 0, 0) 

= 'Vr r)i3 r r))2 (l, 1, 1, 1) 
= r„ )4 r„,3$ 2 (l, 1,1,1) 
= ^,4^,3(1, 1,2,2) 
= r V $s lS2 (l,l,2,2) 
= r^(l,l,2,3) 
= $ Sl s 2 s 3 (l,l,2,3) 
= (1,1,2,4). 

We complete the generation by permuting each component into its correct 
position 

S (J '(?7,2)S( T '(7?,3)S ( j'( ?? ,4)(l, 1, 2, 4) = S<j'( r?) 2)S (J /(?7,3)S3S2(l) 1, 2, 4) 

= S (T '(^ i 2)S (T '( r , > 3)(l,2,4, 1) 

= s CT '(^2)S 2 (l,2,4, 1) 

= 5^(^,2) (1,4, 2,1) 

= si(l,4,2,l) 
= (4,1,2,1). 

We now move onto the major goal of the paper, developing the recur- 
sive generation algorithm for the nonsymmetric and interpolation Macdonald 
polynomials. 

3.3 Recursively generating polynomials 

Beyond the nonsymmetric and interpolation Macdonald polynomials there 
are many other families of polynomials that can be recursively generated via 
switching and raising type operators. Some examples include nonsymmetric 
and interpolation Jack polynomials (see e.g. (TO! Chap. 12]) and the gener- 
alised nonsymmetric Hermite and Laguerre polynomials (see e.g. [TUl Chap. 
13]). In this section we show how the algorithm developed in Proposition 
CD can be employed to recursively generate any of the polynomials in these 
families. 

For simplicity in this section we use F v (z) to denote any of the afore- 
mentioned families of polynomials. We show how the recursive generation 
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algorithm works in the general setting and then provide a specific example 
for the nonsymmetric Macdonald polynomials. 

To most simply express the sequence of operators required to recur- 
sively generate a composition according to Proposition [1] we use the numbers 
1, ... ,n — 1 to represent the allowable switching operators, to represent the 
raising operator and denote the required sequence by R(t/). For Example [T] 
in the previous section we observe that 

R((4, 1, 2, 1)) = {0, 0, 0, 0, 0, 0, 2, 1, 0, 3, 2, 1, 0, 2, 3, 2, 1}. 

Proposition 2 Define 

(11) 

where r](j) represents the composition obtained after j transformations from 
(0, . . . , 0) to rj specified by R(r/) and F SiV (z) and F^ v (z) are obtained from 
F v (z) using known formulas. With initial input RGi(l, (0, . . . , 0), R(t/)i) and 
each subsequent polynomial derived from the previous by entering the newly 
obtained polynomial and composition along with the next number in R(r/) in 
RG the polynomial F v (z) will be obtained after \R(rj)\ steps. 

Proof By Proposition [1] we know that the sequence specified by R{rf) will 
recursively generate the composition r\ from (0, ... ,0). Consequently RGj 
will recursively generate F v (z) from F(o,...,o) {z) ■ □ 

We note that we keep track of the composition labelling the polynomial at 
each stage due to requirements of the formulas transforming F n (z) to F SiV (z) 
and F$ v (z). 

Example 2 To recursively generate a nonsymmetric Macdonald polynomial 
using the methods in the previous proposition we begin by rewriting the re- 
cursive generation formulas for the nonsymmetric Macdonald polynomials 
as 

r^E^z) - ^J-l^ E^z), r/, < rj i+1 

E SiV (z) = { E r,{z), Vi = Vi+i 

(1 ^" (9 ' f))2 (T t E v (z)-p±-E v (z)), Vi > Vi+l! 



0-tS i , v (q,t))(l-t- 1 Si, v {q,t)) \**"V\*J 1-^ („,*)- 

E^{z) =t*^ >1 ^^ q E v {z). 
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To generate the polynomial E(2u(z; q,t) from E^ Qt0 ^(z; q, t) = 1 using RG we 
first compute R{(2, 1)) using Proposition\J\to be R(2, 1) = {0, 0, 0, 1} we then 
proceed recursively 



RG 1 (l,(0,0),0) = { J E (0 ,i)(z),(0,l)} 
RG 2 (E m (z), (0, 1), 0) = {£ ( i,i)(z), (1, 1)} 
RG 3 (£ ( i,i)(20, (1, 1), 0) = {£(i, 2 )(,2), (1, 2)} 
RG 4 (£(i )2 )(z), (1, 2), 1) = {E (2>1) (z), (2, 1)}. 



In the cases where F v (z) is homogeneous, for example the nonsymmet- 
ric Macdonald and Jack polynomials, we can greatly reduce the number of 
operators required to generate polynomials labelled by compositions with no 
zero components by making use of the relationship [16] 



This result allows us to omit the first n x minj?]} zeros from K(f]), forming 
a new set R'(^), and begin our recursive process with 



{(*!... z n )™^\ (minM, . . . , mm{ V }), R'(r])i} 
rather than {1, (0, . . . , 0), R(r?)i}. 

Example 3 Using (Q^j we recursively generate £(2,1) (z; q, t) from (z; q, t) . 
With R'((2, 1)) = {0, 1} we obtain 



3.4 Runtimes and Software 
3.4.1 Algorithm runtimes 

In this section we analyse the performance of the algorithm constructed in 
Proposition [2] when used to generate the nonsymmetric and interpolation 
Macdonald polynomials. 



Fr,+(k n )( Z ) = («!■■■ ZnfF^z). 



(12) 



RG 1 (z 1 z 2 ,(l,l),0) 
RG 2 (£ (1 , 2) (z),(l,2),l) 



{£(!-,)(*), (1,2)} 

{E (2il) (z),(2,l)}. 
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To provide a thorough analysis we must select polynomials of varying de- 
grees of complexity. That is, varying the number of variables, the maximum 
degree and of these polynomials selecting those requiring the least and most 
number of operators to generate. Of the polynomials with n variables and 
maximum degree k those that take the least and most number of operators 
to generate using the algorithm of Proposition [2] are labelled by composi- 
tions (q n ~ r , (q + l) r ) and (k, 0, . . . , 0), where k = qn + r, and require k and 
nk operators to generate, respectively. We note that we have employed ( IT2|) 
into our algorithm for the nonsymmetric Macdonald polynomials and poly- 
nomials labelled by compositions with no zero components will be generated 
more efficiently than what is specified above, for example the nonsymmet- 
ric Macdonald polynomials labelled by rj = (k n ) would be generated almost 
instantly. 

Table [1] shows the runtimes of the generation of different nonsymmetric 
Macdonald polynomials using the recursive generation algorithm of Proposi- 
tion [2j The computer these runtimes were observed on was an iMac 2.4GHz 
Intel core 2 duo processor in version 7.01.0 of Mathematica. 





E v (z) 


KM 


\v\ 


V 


time 


# operators 


time 


# operators 


4 


(0,4) 


0.0308 


7 


0.7792 


7 




(1,3) 


0.0064 


3 


0.0759 


5 




(2,2) 


0.0011 





0.0009 


4 




(3,1) 


0.0177 


4 


0.6081 


6 




(4,0) 


0.0674 


8 


12.3216 


8 




(0,0,4) 


0.48 


10 


4.54 


10 




(1,1,2) 


0.0014 


1 


0.0010 


4 




(2,1,1) 


0.0105 


3 


0.3108 


6 




(4,0,0) 


0.6337 


12 


65.7861 


12 




(1,3,0) 


0.3171 


9 


30.3733 


9 


7 


(0,7) 


0.2677 


13 


34.4630 


13 




(3,4) 


0.0019 


1 


0.001648 


7 




(4,3) 


0.0070 


2 


0.1471 


8 




(7,0) 


0.4744 


14 


68.2642 


14 



Table 1: Runtimes in seconds and number of operators required for the 
recursive generation algorithm 
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We note that it is the inhomogeniety of the interpolation polynomials 
that makes their construction time longer than the nonsymmetric polyno- 
mials. Furthermore, we observe occasions where fewer operators result in 
longer runtimes, for example £(0,0,4) (2) compared with E( 1>3 ^(z). This can 
be explained by the complexity of the polynomials generated in the recursive 
generation process. 
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3.4.2 Mathematica notebook 

We now provide the details of the Mathematica notebook containing the 
polynomials and operators contained within this paper, and more generally 
the author's PhD thesis [3]. The Mathematica notebook, titled SpecialFunc- 
tions.nb, can be found at www.ms.unimelb.edu.au/ wbaratta/index.html and 
ran on Mathematica 7. The original purpose of the Mathematica notebook 
was to efficiently generate nonsymmetric and nonsymmetric interpolation 
Macdonald polynomials to develop an understanding of their known theory. 
Surpassing this motivation it was extensively used throughout the author's 
PhD candidature to assist with conjecture formulation and testing. We now 
present a table containing the key functions defined in the Mathematica 
notebook. The notations used in the table are consistent with those used 
throughout the paper. For those extending beyond the current paper (de- 
noted by a *), we refer the reader to |3J. Note, some functions are not 
introduced until the following section. 



Syntax 


Description 


CompositionModulus [rj\ 


Computes 77 1 


SwitchComposition[?7, i] 


Computes Sifj 


RaiseComposition [rj] 


Computes Qrj 


*c%/] 


Computes ci(r)) 


*Dominance[?7, A] 


Determines whether rj < A or A < i] 


*PartialOrder[?7, A] 


Determines whether rj H A or A ^ rj 


*PartialOrder2[?7, A] 


Determines whether rj < A or A < i] 


*Successor[?7, A] 


Determines whether rj A or A -<' rj 


*ArmLength[?7, 


Computes a v (i,j), the armlength of a composition at a square 


*ArmCoLength[r], i, j] 


Computes a'(i,j), the coarmlength of a composition at a square 


*LegLength[?7, 


Computes l v (i,j), the leglength of a composition at a square 


LegCoLength^, i] 


Computes l' v (i), the coleglength of a composition at a square 


*Md[r]} 


Computes d v (q, t) 


*MdDash[r?] 


Computes d' (q, t) 


*Me[rj] 


Computes e n (q, t) 


*MeDash[r?] 


Computes e' v (q, t) 


Rfo] 


Computes the set R(r/), specifying the sequence of operators 




required to generate rj from (0, . . . , 0) 



Table 2: Mathematica functions relating to compositions 



15 



Syntax 



Description 



PermutationOnComposition[o", 77] 
DecompositionOnComposition [{ , 
SwitchingOperatorf/, i, j] 
PermutationOnPolynomial[(x, /] 
S hort est Permut at ion [77] 

RequiredPermutationfr/, A] 



Computes a(rj) 
J, 77] Computes ... s^r? 

Computes s^/(. . . , ^, . . . , Zj, . . .) 
Computes af(z±, . . . , z n ) 
Computes u v , the shortest 
permutation such that cu~ 1 (t7) = 77 4 
Computes the permutation a such 
that (7(77) = ^ 



Table 3: Mathematica functions relating to permutations 



Syntax 


Description 


* Monomial [77] 


Computes z v 


ElementarySymmetricFunctionfr, n] 


Computes e r (z) in n variables 


*CompleteSymmetricFunction[r, n] 


Computes h r (z) in n variables 


SymmetricMonomialFunction [k] 


Computes m K (z) 


* Vandermonde [n] 


Computes A(z) 


* Vandermonde J [J] 


Computes A J (z) 


t Vandermonde [n] 


Computes A t (z) 


*t Vandermonde J [J] 


Computes Af(z) 


Schur[/t] 


Computes s K (z) 


Zonal [k] 


Computes Z K (z) 


HallLittlewood[/t] 


Computes P K (z; t) 



Table 4: Mathematica functions relating to miscellaneous polynomials 
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Syntax 


Description 


SymJack[/t] 


Computes P K (z; a) 


NSJack^] 


Computes E v (z; a) 


IntJack[r/] 


Computes E*(z;a) 


*JEvalue[r?] 


Computes rj a , where Tj i := arji — l' n (i) 


*EOpSymJack[/, n\ 


Computes D 2 (ot)f 


*EOpNSJack[/, n, i] 


Computes 


*EOpInt Jack[/, n, i] 


Computes Sf/ 


*JInnerProduct[/, g, n, k] 


Computes (f,g)i/k for polynomials /, g of n variables 


Table 5: Mathematica functions relating to Jack polynomials 


Syntax 


Description 


NSMacfa] 


Computes E v (z; q, t) 


SymMac[re] 


Computes P K (z; q, t) 


ASymMac[/t] 


Computes S K (z; q, t) 


IntMacfa] 


Computes E*(z; q, t) 


SymIntMac[fi;] 


Computes P* (z; q, t) 


ASymlntMac [k] 


Computes S*(z; q, t) 


MEvalue[r?] 


Computes fj, where r] i := q^H 1 '^ 


*EOpSymMac[/, n\ 


Computes D\(q,t)f 


*EOpNSMac[/,i] 


Computes Yif 


*EOpIntMac[/, i] 


Computes Sj/ 




Computes T.J 


Tilnv[/, i] 


Computes T~ x j 


Phiq[/, n] 


Computes $ 9 / 


Hi[/,i] 


Computes Hif 


Phiqlnt [/, n] 


Computes $*/ 



*MInnerProduct[/, g, n, k] Computes {f,g) q q k, for polynomials /, g of n variables 



Table 6: Mathematica functions relating to Macdonald polynomials 
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Syntax Description 



*UPlus[/, n] Computes U + f, where / is a function of n variables 

*UMinus[/, n] Computes U~f, where / is a function of n variables 

*UPlusInt[/, n] Computes U+f, where / is a function of n variables 

*UMinusInt[/, n] Computes U~f, where / is a function of n variables 

OIJ[/, n, I, J] Computes Oi t jf, where / is a function of n variables 

OIJInt[/, n, I, J] Computes Oj^f, where / is a function of n variables 

PreSymMacfrf , T, J] Computes Sfi J \z; q, t) 

PreSymIntMac[r]*, /, J] Computes S*l (I ' J \z; q, t) 

aEtafr], T, J] Computes the normalisation alf' J \q,t) 

Table 7: Mathematica functions relating to prescribed symmetry polynomials 



Syntax 


Description 


*JEta[7?,r] 


Produces the set of A such that A = \r] \ + r 


and 


V =<' A =<' v+ (i n ) 


Pieri[?7, A, r] 


Computes the Pieri-type coefficient A^{(q,t) 


*GeneralisedBinomialCoemcient[A, rf\ 


Computes (^) ^ 



Table 8: Mathematica functions relating to Pieri-type coefficients 



Copyleft Copyleft 2011 Wendy Baratta 

Permission is granted to anyone to use, modify and redistribute SpecialFunction.nb 
freely subject to the following. 

• We make no guarantees that the software is free of defects. 

• We accept no responsibility for the consequences of using this software. 

• All explicit use of this notebook must be explicitly represented. 

• No form of this software may be included or redistributed in a library to be sold 
for profit without our consent. 
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4 Further Work 



It is hoped that the notebook SpecialFunctions.nb will continue to assist 
researchers with conjecture formulation and testing in areas relating to the 
polynomial families discussed within this paper. In fact this has already 
shown itself to be the case in a study of special vanishing properties of Jack 
polynomials for a = — (r — l)/(k + 1) and Macdonald polynomials with 
t h+l q r ~ 1 = 1 [5]. Below are some examples of computations that aim to 
provide motivation for two research problems stemming from the author's 
previous studies; Macdonald polynomials with prescribed symmetry [2] and 
Pieri-type formulas for nonsymmetric Macdonald polynomials [3]. 

4.1 Interpolation polynomials with prescribed symme- 
try 

In [2] the author investigated properties of Macdonald polynomials with pre- 
scribed symmetry. Macdonald polynomials with prescribed symmetry, de- 
noted S^* J \z\ q,t), are generalisations of the symmetric Macdonald polyno- 
mials and the antisymmetric Macdonald polynomials, denoted S v + + s(z; q, t) 
where 5 = (n — 1, . . . , 1, 0), and are symmetric with respect to some vari- 
ables and antisymmetric with respect to others. They are obtained from 
the nonsymmetric Macdonald polynomials via symmetrisation by Oi t j and 
normalisation. The subscript /, J indicates the sets of variables for which 
Oij symmetrises and antisymmetrises with respect to. The operator Oj j is 
defined by 

where l(u) is the length of the permutation u, 

T u := T h . . . T iv (13) 

where . . . is a reduced decomposition of u, and Wi u j := (sk, k G IU J), 
a subset of S n where each u G Wj u j can be decomposed as 

us = UjUj, with cu j G Wj and Uj G Wj. 

The prescribed symmetric Macdonald polynomials is obtained from the 
nonsymmetric Macdonald polynomials according to 

O u E ri {z-q ) t) = a^S^ J \z-q ) t) 1 
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where if is a composition satisfying 

Vi > Vi+i f° r alH G / and rj* > r]* +1 for all j G J 

and alj'^ is a normalisation that ensures the coefficient of z 71 * := z^ 1 . . . zH n 
is unity. 

The theory of Macdonald polynomials with prescribed symmetry leads 
most naturally to an investigation of interpolation polynomials with pre- 
scribed symmetry. Since the operator Hi (J3J) satisfies the Hecke algebra (J2J) 
and plays the same role in interpolation theory as Tj does in nonsymmetric 
theory it seems natural to define the prescribed symmetry operator for the 
interpolation polynomials as 

( " t) 

uiGW iuj 

where H u is as in (fl3l) . 

Trial computation suggest that it may be possible to extend many results 
obtained in [2] to their interpolation polynomial analogues (see for example, 
Propositions 3.2.1, 3.2.2 and 3.3.1). 

Specifically, Proposition 3.3.1 in [2] specifies the relationship between the 
antisymmetric and symmetric Macdonald polynomial 

S v++S {z; q, t) = A t (z)P v +(z;q,qt), 

where A t (z) is the t-Vandermonde and specified by 

A t (z):= J] (^-rV). 

l<i<j<n 

This property has extensions in the theory of Macdonald polynomials with 
prescribed symmetry and is therefore a suitable starting point for the theory 
of interpolation Macdonald polynomials with prescribed symmetry. 

Here we provide some explicit computations obtained using the Mathe- 
matica notebook that may assist with the identification of the relationship 
between the antisymmetric and symmetric interpolation polynomial. 
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In[l] SymIntMac[{l,0}] 

t (X\ + X 2 ) — t — 1 

Out [1] -i-i ^ 

t 

In [2] ASymIntMac[{2,0}] 

(xi - tx 2 ) (t (x\ + x 2 ) —qt—1) 



Out [2] 



Out [2] 



t 



In[l] SymIntMac[{l, 1}] 

r _ (tX\ — 1) (tXo — 1) 

Out [1] ^ ^-^ i 

In [2] ASymlntMac [{2, 1}] 

(xi - tx 2 ) {txi - 1) (tx 2 - I) 



t 2 



In[l] SymlntMac [{1,0,0}] 

t 2 (xi + x 2 + x 3 ) - t 2 - t - 1 
Out [1] - 

t z 

In [2] ASymlntMac [{3, 1,0}] 

Qut r 2 -, (xi - tx 2 ) Ogi ~ fafr) O2 - tx 3 ) (t 2 (xi +x 2 + x 3 ) - qt{qt + !)-!) 

t 2 

These computations show a clear relationship between S* +S (z) and P^(z), 
highlighting that the difficulty in identifying the relationship lies in the trans- 
formations of the parameters. It also appears that the t-Vandermonde prod- 
uct that relates the polynomials is not of standard form. 

4.2 Pieri-type formulas for nonsymmetric Macdonald 
polynomials 

In [3] the formulas for the coefficients in the expansion of a nonsymmetric 
Macdonald polynomial and an elementary symmetric function were given. 
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Explicitly, the formulas for in 

e r (z)E v (z;q-\r 1 )= £ A^E\(z\ q" 1 , t -1 ), 

A:|A| = |r;|+r 

where 

e r (z) := ^ ] • • • Zi n - 

l<il<. ,,<i T <n 

The formulas of A^ given in [3] are given in summation form, though compu- 
tational evidence suggests that almost all of the coefficients display a product 
structure. Here we present explicit formulas for the Pieri-type coefficients for 
the case r = 2 showing which cases display a product structure and which 
cases do not. 



In[l] Pieri[{l,0,l,0},{2,0,2,0},2] In [2] Pieri [{1, 0, 1, 0}, {2, 1, 1, 0}, 2] 

q-1 



Out [1] 1 Out [2] 



qt- 1 



In [3] Pieri [{1,0, 1,0}, {1,1, 2,0}, 2] In [4] Pieri [{1, 0, 1, 0}, {2, 0, 1, 1}, 2] 

out [3 ] '(g-Dfr 3 ,- 1 ) out M ifa-iH* 3 -!) 

(qt 2 - I) 2 (qt 2 - l) 2 



In [5] Pieri [{1,0, 1,0}, {1,1, 1,1}, 2] In [6] Pieri [{1, 0, 1, 0}, {1, 2, 0, 1}, 2] 

Out [5] (g-l)(*-P(f-l) 0ut[6] g^fa -D(*-l)' 

(qP-lf (qt - 1)* - I) 2 



In [7] Pieri [{1,0, 1,0}, {1,2, 1,0}, 2] In [8] Pieri [{1, 0, 1, 0}, {2, 1, 0, 1}, 2] 

Out [7] -^M 0ut[B] ^-D(*-D a 



(qt-lY(qt*-l) (qt-iy(qt*-l) 
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In [9] Pieri[{l,0,l,0},{0,2,2,0},2] In[10] Pieri [{1, 0, 1, 0}, {2, 0, 0, 2}, 2] 

qt 2 (q-l)(t-l) qt 2 (q-l)(t-l) 
Out [9] t --7 --7 — = — Out [10] 



(qt - l)(qt + 1) (qt 2 - 1) (qt - l)(qt + 1) (qt 2 - 1) 



In [11] Pieri [{1,0, 1,0}, {1,0, 2,1}, 2] 

(g - l)t 2 (g (t (g ((g - 2)gt 5 + (3g - l)t 4 + (2 - 3g)t 3 - 2t + 1) + 3(t - 1)) + 2) - 1) 

Out [11] 

(gt - 1) {qt 2 - l) 2 (g 2 i 3 - 1) 

(2) 

In the above formulas for the Pieri-type coefficients A^ 1Q10 ^ x (q,t) a simple 
product structure in q and t, with roots in t being simple fractional powers 
of q for example, is exhibited in all cases except A = (1,0,2, 1). Additional 
trial computations suggest that it is always the successor of the form $$77 
that cannot be expressed simple product. 
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